Improving quantum interferometry by using entanglement 
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! Abstract 
I . 

I ' 

^ ' We address the use of entanglement to improve the precision of generahzed 

' quantum interferometry, i.e. of binary measurements aimed to determine 

whether or not a perturbation has been appUed by a given device. For the 
most relevant operations in quantum optics, we evaluate the optimal detection 
strategy and the ultimate bounds to the minimum detectable perturbation. 
Our results indicate that entanglement-assisted strategies improve the dis- 
crimination in comparison with conventional schemes. A concrete setup to 
approach performances of the optimal strategies is also suggested. 

^ I. INTRODUCTION 

Any interferometric setup is devised to reveal minute perturbations to a given configu- 
^ . ration. Such perturbations may be induced by the environment or by the action of a given 
device. In an interferometer, the internal quantum operation is monitored by probing the 
output state, which, in turn, results from the evolution of a given input. By suitably choos- 
ing the input signal and the detection stage one optimizes the interferometric measurement. 
Optimization has two main goals: i) to maximize the probability of revealing a perturba- 
tion, when it occurs, and ii) to minimize the value of the smallest perturbation that can be 
effectively detected. 

In essence, any interferometric scheme may be viewed as a binary communication system 
mj^, with the perturbation playing the role of the encoded information. In order to see better 
this analogy let us consider the scheme shown in Fig. |l|a. A source S of quantum state 
prepares the input signal, say qq, which travels along the interferometer, and it is eventually 
measured by some detector, denoted by D. The detector is described by an operator-valued 
probability measure (POVM) n(x), with x G X, X being the manifold describing the 
possible detection outcomes. Inside the interferometer we a have generic quantum device, 
which may or may not perturb the signal, i.e. it performs the quantum operation described 
by the unitary U\. If a perturbation occurs the signal is modified and, at the output, we 
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have the state gx = Ux QqUI- The aim of the detection stage is to discriminate between go 
and its perturbed version qx- An optimized interferometer is a device that is able to tell 
which g, for A as small as possible. Posed in this way, interferometry is naturally viewed as 
a binary decision problem, and the detection stage can be described by a two-value POVM 
{Ho, IIa = I — Ho}, which corresponds to the two possible inferences. 

The main goal of the present paper is to demonstrate the benefits of entanglement in 
binary interferometry. We will show that distinguishability of the two hypothesis (Hq: 
nothing happened and Tix'- a perturbation has occurred) can be improved by: i) using an 
input signal which is entangled with another subsystem, and ii) measuring the two systems 
jointly at the output of the interferometer (see Fig. |l]b). 

In order to optimize the detection strategies, and to show the benefits of entanglement, 
we will make use of results and methods from quantum detection theory applied to binary 
decision This approach is particularly useful for our purposes, since it does not refer 

to any specific detection scheme for the final stage of the interferometer, but rather, owing 
to its generality, it allows to find the ultimate quantum limits to interferometry for specific 
classes of quantum signals. 

In the next Section, in order to establish notation, we briefly review the Neyman-Pearson 
approach to quantum binary decision, and state a lemma about minimum input-output 
overlap. Then, in Section |T| we apply these results to the interferometric detection of per- 
turbations induced by the most relevant operations in quantum optics such as displacement, 
squeezing, mixing and phase-shifting. As we will see, entanglement-assisted interferometers 
provide better discrimination than conventional schemes. In Section ^ we analyze an in- 
terferometric configuration that achieves, for the quantum operations discussed in Section 
lT|, the ultimate bounds to precision. Finally, in Section we close the paper with some 
concluding remarks 



II. QUANTUM BINARY DECISIONS IN THE NEYMAN-PEARSON 

APPROACH 

The problem that we are facing is to decide among two hypotheses Hq and Hi about 
the state of a system, which is described by a density operator g on the Hilbert space of the 
system. To each hypotheses it will correspond a different density operator as follows 

Ho '■ the system is in the state go, 
Hx '■ the system is in the state ^a- 

Of course, there are many different measurements which can provide information about the 
state of the system: each of them, however, can be recast mathematically as a two-value 
POVM, corresponding to the two possible inferences Hq and Hi, namely 

Ho, Ha > Ho + Ha = I. (2) 

One then needs an optimization strategy in order to determine the most reliable measure- 
ment discriminating between the two states. If go and gx are orthogonal i.e. gogx = QxQo = 
the solution is trivial, since Ho is the projection into any subspace that contains the support 
of go and is orthogonal to the support of gx, and IIa is simply the complement IIa = I — IIq. 
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In most cases of interest, however, the states are not orthogonal and one has to apply an 
optimization strategy. Since interferometric schemes are frequently used for detecting low- 
rate events, we may want to look for a strategy that keeps a low-rate of false alarm, namely 
of wrong inference of perturbation occurrence. For this purpose, it is suitable to adopt a 
so-called Neyman-Pearson (NP) detection strategy, which consists in fixing a tolerable value 
of the false-alarm probability Qq — the probability of inferring that the state of the system 
is while it is actually go — and then maximizing the detection probability Qx, i.e. the 
probability of a correct inference of hypothesis Tix 0] . It has been proved by Helstrom 
and Holevo [Q that this problem can be solved by diagonalizing the operator 

Qx - f^Qo , (3) 

/i playing the role of a Lagrange multiplier accounting for the bound of fixed false alarm 
probability. According to the optimal POVM is the one in which 11^ is the projection onto 
the eigenspaces of (|]) relative to positive eigenvalues and Hq = I — 11^. Unfortunately, the 
diagonalization of @ is generally not easy. However, when qq = \ipo){i^o\ and gx = \ipx){i^\\ 
are pure states it can be easily solved analytically, by expanding |?/'o) and \ipx) on the 
eigenvectors of the difference operator (|^). In this way one can evaluate both Qq and Qx 
versus /i, and after eliminating /i from their expressions one obtains 



Q> 



for < Qo < l/tp, 



v/Q^+v/(l-Qo)(l-|«:P)^ 

for |Kp < Qo<l. 

where |kP = | (t/'oIV'a) P = I (V^oI^^aI'^/'o) P is the overlap between the two states. The detection 
probability is a decreasing function of the overlap — the smaller the overlap, the easier the 
discrimination — since one can reach detection probability 1 while keeping a low false alarm 
probability. On the contrary, when the overlap approaches 1 one is forced to decrease the 
detection probability in order to keep the false alarm probability small. 

The optimal choice of the probe that minimizes the overlap, depends on the eigenvalues 
of the unitary operation Ux- In order to illustrate this, let us expand Ux in terms of its 
eigenvectors Ux = J2j ^^"^^ (with integrals replacing sums in case of continuous 

spectrum) and let's denote by 0{Ux) = min^ | ('?/'|?7a|V') P the minimum overlap between the 
two possible outputs, as obtained by varying the probe state. Then we have the following 
overlap Lemma 0,0: the minimum overlap 0{Ux) is given by the distance from the origin 
in the complex plane of the polygon whose vertices are the eigenvalues of Ux- Therefore, the 
overlap is either zero (if the polygon includes the origin) or it is given by 

0{Ux) = cos' ^ , (5) 

where A(p is the angular spread of the eigenvalues. Zero overlap can be achieved with a probe 
state that is given by a superposition of at least three eigenvectors of Ux, corresponding to 
eigenvalues that make a polygon that encloses the origin (or, if they exist, by a superposition 
of two of them corresponding to diametrically opposed eigenvalues). Instead, if the minimum 
overlap is not zero, it is achieved by the optimal probe state given by 

\^) = ^{\^^) + \^,)) , (6) 



with A(p = ipi — (pj. 
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III. ENTANGLEMENT IN QUANTUM INTERFEROMETRY 



In this section we compare the performances of single-mode (Fig. |I]a) and entanglement- 
assisted interferometric schemes (Fig. |l|b) in the detection of small perturbations induced 
by relevant quantum optical operations such as displacement, squeezing, mixing and phase- 
shifting. The comparison is made in terms of the detection sensitivity, namely, upon 
parametrizing the "size" of the perturbation — whence the corresponding output state — by a 
coupling parameter A. In other words, the comparison is made in terms of the minimum de- 
tectable value Amin of A corresponding to output states that can be effectively discriminated 
while keeping the acceptance ratio 7* of the NP strategy large, namely 7* = Qx/Qo ^ 1. 
We will call the quantity Xmin the "sensitivity" of the interferometric scheme. Using Eq. (^ 
the above condition can be written in term of the overlap as follows 

|/t|2 = l-A(Qo,7'^) (7) 



A(Qo, 7*) = Qo[l+ I'^il - 2Qo) - 2 V7'^(l - Qo){l - 7*Qo) 

For each class of transformations, we will make some general considerations and then focus 
our attention on sensitivity bounds that can be achieved using realistic {i.e. feasible with 
current technology) probe signals. 



A. Perturbation made of a single-mode complex displacement 

Let us first consider the case when the perturbation is imposed by the displacement 
operator Ua = D{a) = exp(aa^ — aa). In principle, in this case, the discrimination can 
be done exactly with single-mode probe. This can be seen by writing the displacement as 
Ua = exp{i2\a\xe), xe = l/2(a^e*^ + ae~'^) being the quadrature operator, and 9 = arg(Q!). 
Since the spectrum of the quadrature coincides with the real axis, the spectrum of Ua 
covers the whole unit circle, and, therefore, the states {ipo) and {ipa) = f^al^o) can be 
discriminated with certainty either by choosing I^/'q) as the eigenstate of the conjugated 
quadrature xe+T,/2, or, according to the overlap lemma, as a superposition of at least two 
eigenstates of the quadrature xq. Unfortunately, such optimal states are unphysical, since 
they are not normalizable and have infinite energy. Moreover, even though we approximate 
them with physical states with finite energy, the identification of the optimal states would 
require the knowledge of the phase of the perturbation. In order to see that, let us rewrite 
the eigenvector 1 0)0+^/2 as the limiting case of a squeezed vacuum, |0)e)+,r/2 = lim|(^|^oo |C) = 
lim|^l^oo 'S'(C)lO), where 9 = arg(C) -|- 7r/2 is the argument of the squeezing parameter ( of 
the squeezing operator given by 

5(C) = exp[l/2(C^a+2-CV)], (8) 

and |0) is the electromagnetic vacuum. Our squeezed vacuum has mean photon number 
N = sinh^ \(\. The overlap is readily evaluated as 

IkI^ = |(C|£'(a)|C)r = exp|-|a|2 2N + 1 + ^/ N{N + 1) cos26 } , (9) 
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where S = arg((^) — arg(a). By inserting the overlap in Eq. (|^) we obtain the minimum 

on the phase 
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detectable jap. However, Eq. (P) shows a very strong dependence of 
parameter 6, which makes the whole optimized scheme very unstable, namely one should 
know the phase of perturbation very precisely in order to get a truly optimized detection. 
Indeed, we have 



\a\ 



mm 
2 A^>1 



4NA{Qo,Y 



for 6 
for 6 



7l/2 

, 



(10) 

(11) 



with the second expression that shows an asymptotically divergent behavior in A^. 

Let us now consider an entanglement-assisted scheme. We suppose you have available 
a two-mode probe state and consider the configuration Ua = D{a) ® I in which the 
displacement perturbs one mode, say a, and the other mode is left unperturbed. As the 
probe state we consider the entangled state from parametric downconversion of vacuum for 
finite gain — the so-called "twin-beam" state 



= Vl - ^ x" \nn)) < X < 1 



(12) 



where here \nn)) = \n)a (S> \n)b. The twin-beam in Eq. ([T^) has mean photon number 
A^ = 2x^/(1 — x^) and it is achieved starting from the vacuum via the unitary evolution 
|x)) = exp[x(a^6^ — ab)]\0)). In order to evaluate the sensitivity, the main task is now to 



calculate the overlap \k\ 



(x\Un\x) 



We have 



X 



x™+" {{mm\D{a) ® I\nn)) 



m=0 n=0 
oo 



:i-x' 



exp 



\a 



J2 a;'" {n\D{a)\n) = (1 - x') e"^!"!' 

n=0 
2l + x2 



n=0 



2 l-x2 
Equation (|T3|) implies for \a 



exp 



\a\ 



■{N+l) 



(13) 



2 

min 



the scaling 



\a\ 



N+l 



(14) 



which is remarkably independent on the phase of perturbation, and thus represents a robust 
bound to the sensitivity of a single- mode displacement. 



B. Perturbation made of a single-mode squeezing (phase-sensitive amplifier) 

The second kind of perturbation that we analyze is the squeezing of a single radiation 
mode, which is described by the squeezing operator S{() in Eq. (|). Without loss of 
generality we can consider ( = ( = r a.s real and use the notation Ur to indicate the 
transformation, namely Ur = exp[—irA], with A = z/2(a+^ — a^). The spectrum of A is 
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continuous and extends from — oo to oo: this means that the eigenvalues of Ur cover 
the whole unit circle. Therefore, it is possible in principle to discriminate the perturbation 
exactly, using as a probe either an eigenstate of the operator conjugated to A, or using 
a superposition of two or more eigenstates of A. However, analogously to the case of the 
displacement, such probe states would be non normalizable and have infinite energy, whence 
one must resort to physical approximations of such states. For a coherent probe the overlap 
can be calculated through the overlap of the corresponding Wigner functions, giving as a 
result 



(a|f/,.|a)p = 

2N cos^ — cosh r — sinh r)^ 2N sin^ 0(1 — cosh r + sinh r) 

exp 



1 + exp(2r) 1 + exp(— 2r) 



(15) 



where = jap = (a|a^a|a) is the mean number of photons of the probe state. By expanding 
for small r we have 

|(a|f/,|a)P~ l-A^r^, (16) 
and therefore the minimum detectable perturbation would be 



A(Qo,r) 



(17) 



For a squeezed vacuum probe S'(C)|0) one has 

K= (0|5l'(C)?7r5(C)|0) = [cosh r + 2i sinh Id cosh Id sinh r sin V-]-^ (18) 
where ip = arg(d and correspondingly the minimum detectable r is given by: 



log[T3A(kr) 1 - VA(Qo,r)(2-A(go,7*)) 



for sinip = 0, 



^ , otherwise, 

2 N sin If) ' 



(19) 



with N = sinh^ |d- The bound in Eq. ([T9| ) strongly depends on the phase between the 
squeezing perturbation and the squeezing of the probe, and therefore cannot be achieved in 
practice without prior knowledge of the phase of the perturbation. 

Let us now consider an entangled probe state in a twin beam state of the form (^). The 
input-output overlap is calculated as follows 



K 



{{x\Ur ® I|a;)) = (1 - x^) J2 x^"{n\Ur\n) . (20) 



n=0 



In order to calculate the matrix element {n\Ur\n) we use the identities S{r) 

tanh(r)at |^cosh(r)] ^e~ ^ tanh(r)a g^-^d 



e-|tanh«a|^^ = ^^^7^1^ " 2/), (21) 
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where [m] indicates the integer part of m, and finally we get 



/ irr , V n\ ^ (-l)nsinh(r)]2' 

(n\Ur\n) = ^> ^tt-Ut ^^rr- (22) 



Using Eq. (^21) we calculate k by means of Eq. (|2 



CX) 



LCOsh(r)r+^ 4H/!)2(n-20! 

~ [cosh(r)]i V 4cosh2(r); F ^ n!2/! 

[(x^ + l)cosh(r) - 2a;2]5' 
Inserting this expression in Eq. (|^) we have for Tmin the scaling law 



^' l-A(Qo,7^) ViV2 + 2iV + 2 V l-AWo,7^)iV ^ ' 



The same result is obtained by varying the phase of the squeezing amplitude C, thus con- 
firming the robustness of the bound (p3D that is obtained using an entangled probe. 



C. Perturbation made of a two-mode phase-shift 

The third problem we address is that of a perturbation induced by the two-modes phase 
shift operator a^h + ah\ characterizing a mixer (beam splitter) or a Mach-Zehnder inter- 
ferometer. This case differs from the previous ones in that the perturbation is represented 
by the two-modes unitary operator = exp{i0(a^6 + ab^)}. In this case the spectrum is 
given by exp{im0}, with m G Z (see e.g. [0). Therefore, if = {q/p)7i with q E 2Z + 1 
and p G Z (but this is a null-measure set of values of 0) then the optimal state is given 
by a superposition of two eigenstates of with eigenvalues differing by vr. In the general 
case, the optimal state is any superposition of three or more eigenstates of V^, such that the 
polygon of its eigenvalues on the unit circle encloses the origin 0. Such optimal states are 
entangled, since they are obtained from the eigenstates |n, d)) of a^a — b'^b 

(au-b%M)=,M), in,rf»={j:);''|i;'i|) (24) 

by the unitary transformation exp{ — (7r/4)(a^6 — a6^)}. Actually, the optimal states are far 
from being practically realizable. However, we have proved that they are entangled, and 
this suggests to explore the possibility of performing a reliable discrimination by physically 
realizable entangled states. For a twin-beam we have 

K = {{x\V^\x)) = 

= (1 - a;2)((00|e^'''V^°"^''e^57i{«ta-fet6)}gi{T,oa6t}gX.at6t|QQ^^^ ^25) 
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where 70 = tan and 71 = 



— log(cos^0). After some algebra we get 



\K 



2 



1 + Y^ifJ 1 + N{N + 2) sin^ ■ ^^^^ 



The minimum detectable 0, according to (p6D, is thus given by 



. - arcsin ( J^^iXL] - ^^^"'^^^ (27) 
0_-arcsm|^^____J _ ^ . (27j 

The scaling in Eq. (pTj) does not depend on any parameter but the energy of the input state. 
This should be compared with the sensitivity of the customary single-mode interferometry 



ri| based on squeezed states, where the same scaling is achieved only for a very precise tun- 
ing of the phase of the squeezing. This means that the entanglement-assisted interferometry 
provides a much more reliable and easily tunable scheme. 



IV. IMPLEMENTATION BY DIFFERENCE-PHOTOCURRENT 

INTERFEROMETRY 

In this section we present a concrete scheme for binary decision based on an entangled 
probe. The scheme should be feasible with current technology, and would allow to approach 
the ultimate precision bounds that have been obtained in the previous sections. In Fig. |^ 
we show a schematic diagram of the interferometric setup. The input state is the entangled 
twin-beam \x)) produced by a nondegenerate optical parametric amplifier (NOPA). Such 
entangled probe is possibly subjected to the action of the unitary U\ (Figs. |]a and ^jb 
describe the cases of a single mode and a two- mode perturbation respectively). At the output 
the two beams are detected and the difference photocurrent D = a^a — b'^b is measured. If 
no perturbation occurs, then the output state is still a twin-beam, and since \x)) is an 
eigentstate of D with zero eigenvalue we have a constant zero outcome for the difference 
photocurrent. On the other hand, when a perturbation occurs the output state is no longer 
an eigenstate of D, and we detect fluctuations which signal the presence of the perturbation 
itself. The false-alarm and the detection probabilities are given by 

Qo = p(d ^ O|not Ux) = (28) 
Qx = P{d ^ 0\Ux) = 1-P{d = 0\Ux) , (29) 

where the probability of observing zero counts at the output, after the action of Ux, is given 
by 

P{d = 0\Ux) = Y, K(^'^I^a|x))|' , (30) 

n 

since the eigenvalue ci = is degenerate. In this case the false-alarm probability is zero and 
therefore it is not necessary to introduce an acceptance ratio. The scaling of the minimum 
detectable perturbation can be obtained directly in term of the detection probability Qx by 
Eqs. (|29|) and (0). For the three transformations considered in the previous Section we 
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have 

Displacement 

P(d = 0\a ^ 0) = exp + A^)) /o(|a| V^(^ + 2)) 



\a\ 



N 



P{d = 0\r ^ 0) 



P(c/ = 010^0) 



Squeezing 

+ 0{r^) — > 1 
Two-mode phase-shift 



N 



V2Q> 
N 



(31) 



(32) 



(33) 



On can see that in all examples considered above, a realistic interferometer based on a 
difference-photo current measurement provides a precision that re-scales with the energy in 
the same way as the ultimate bounds obtained in the previous sections. 

It is worth noticing that the experimental measurement of a modulated absorption based 
on entanglement-assisted difference-photocurrent detection has been already performed us- 
ing the entangled beam exiting an amplifier above threshold (parametric oscillator, OPO) 

Hi- 



V. SUMMARY AND CONCLUSIONS 

In this paper we have analyzed the effect of entanglement on the interferometric estima- 
tion of relevant quantum optical parameters such as displacing and squeezing amplitudes 
or interferometric phase-shift. We have evaluated the minimum detectable perturbation 
according to the Neyman-Pearson detection strategy, and have shown that entanglement 
always improves the detection in comparison with conventional schemes. In particular, for 
the case of estimation of the displacement and the squeezing amplitudes we have shown 
that the precision of the apparatus that use an entangled probe is independent of the phase 
of the perturbation, and is therefore more robust and reliable of a simple scheme based on 
single-mode probe states. Similarly, for the estimation of a two-mode phase-shift, entangle- 
ment the interferometer is much more stable when we use a twin beam that when we use 
squeezed states. 

Since the Neyman-Pearson detection strategy does not correspond to a realistic detec- 
tor, we proposed a feasible interferometric setup that is based on the measurement of the 
difference photocurrent on an entangled twin-beam. Remarkably, this scheme has the same 
energy-scaling of the ultimate precision bound, and at the same time is very stable. We 
conclude that the technology of entanglement can be of great help in improving precision 
and stability of quantum interferometers. 
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FIGURES 




FIG. 1. A generalized interferometer is a binary detection scheme aimed to check whether or 
not a given quantum device (the hexagon in the figure) has performed the quantum operations 
described by the unitary operator Ux. The signal employed as a probe is prepared by the source 
S and then enters the device, which may or may not apply Ux. The two hypothesis: TCq (the 
signal is unperturbed) and 7ix (Ux has been applied) should be discriminated on the basis of the 
outcome of the detector D. (a): simple scheme involving a single-mode probe, (b): scheme involving 
entanglement-assisted binary detection. 




(b) 

FIG. 2. Interferometric scheme to achieve ultimate bounds on precision by means of an entan- 
gled probe. The NOPA generates a twin-beam which may be subjected to the action of the unitary 
Ux- At the output the beams are detected and the difference photocurrent is measured. For an 
unperturbed interferometer the output is again a twin-beam state, and the scheme is designed in 
order to obtain a constant zero difference photocurrent, whereas a perturbation Ux would produce 
fluctuations in the difference photocurrent. (a) general scheme for single-mode perturbation; (b) 
scheme for two- mode perturbation. 
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